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the controller variable that is coupled to the controlled system. This causes decoherence in the 
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equations are similar in structure to recently proposed models for consistently combining quantum 
and classical stochastic dynamics. 
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I. INTRODUCTION 

The results of measurements made on quantum systems are classical random variables extracted from classical 
instruments. A great deal of effort has been spent in trying to reconcile this fact with the apparent universality 
of quantum theory. In one approach, referred to as the decoherence program[I], the open, dissipative nature of a 
measurement apparatus becomes central. An equivalent, although often ignored, problem of consistency arises when 
we ask for the conditions under which the control parameters in a quantum system can be regarded as entirely 
classical: why not simply include the control apparatus as part of the quantum system as well? In this paper it will 
be shown that the classical controllability of a quantum system also requires us to make explicit reference to the open 
and dissipative nature of control devices (lasers for example). 

The question of specifying the conditions under which we can treat the control of quantum systems using classi- 
cal signals has arisen in a variety of fields. Most recently it has engendered some debate in the field of quantum 
information [2j [3]. Quantum computing requires one to be able to implement rather complex time dependent Hamil- 
tonians for many interacting quantum system( qubits). If any of these subsystems became entangled with the control 
apparatus, it would not be possible to perform unitary control and errors would necessarily enter the quantum compu- 
tation. For example, in ion trap quantum computing we do not concern ourselves with fundamental quantum nature 
of the controlling laser fields, or the RF trapping fields which are treated entirely classically. A very elegant discussion 
of this problem has been given by Carmichael and co workers[H[S] who draw attention to the fact that lasers are open 
quantum systems of a very special kind and the conditions for classical controllability turn out to be very similar to 
the conditions for the classical stochastic nature of measurement results. In the emerging field of hybrid engineered 
quantum systems we will need to consider a greater variety of quantum controllers than simply a laser. 

There is another control problem in which the control is intrinsically classical and cannot be treated as an appropriate 
limit of a quantum interaction; when the control is implemented by a gravitational field. This is related to the long 
debated inconsistencies that arise when we combine classical and quantum dynamics [6 . One resolution has been 
proposed by Diosi[7]. It requires that we add noise to both the classical and the quantum dynamical system so that 
the noise spectral density in both cases are linked. Adding noise to a classical dynamical system can be quite harmless, 
but adding noise to a quantum system necessarily implies decoherence. 

Wiseman and Warszawski 8 discuss a related measurement problem that requires one to combine classical and 
quantum stochastic dynamics. They consider realistic measurement scenarios in quantum optics in which the actual 
observed stochastic process is conditioned not only on the quantum nature of the source but also by additional classical 
noise processes occurring in the detection circuits. The resulting equations are similar in structure to those obtained 
in this paper, as well as those postulated to provide a consistent hybrid classical and quantum dynamics. A similar 
motivation was used by Reginatto and Hall [5] to consistently combine classical and quantum stochastic processes. 
An interesting application of their model was recently given by Chua et al. [TU] . A different kind of classical-quantum 
hybrid dynamics is discussed by TsangfTT] where the objective is to optimally estimate the classical, and possibly 
stochastic, controls on a quantum system by making continuous measurements on the quantum system. The starting 
point of this paper is different from these models. Here we seek the conditions under which two interacting quantum 
systems can be regarded not as a measurement scenario but rather as a classical control scenario. 

A simple example will suffice to illustrate the relationship between classical controllability and the role of decoherence 
in measurement. Consider the case of applying a classical impulsive force to a free particle, that is to say, a force 
applied on a time scale much shorter than the time scale of the intrinsic dynamics of the particle. The change in the 
state of the system due to the impulsive force is well described by the unitary transformation |Vw) = U(X)\tp in ), 
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where 

U(X) = e~ iX ^ n (1) 

where X is a classical control parameter equal to the change in momentum of the particle due to the impulse and 
where (q,p) are the quantum canonical position and momentum operators of the quantum system. 

Now suppose the impulse was in fact due to a collision with another free particle — call it the control particle for 
ease of reference — with canonical position and momentum operators (Q,P). In this case we know there will be a 
mutual exchange of momentum and the change in the state through the collision may be described by the unitary 
operator 

U = e r lK ^ /h (2) 

Under what conditions are we justified in describing this interaction by the unitary transformation in Eq.Q? 

The answer of course depends on how we prepare the state of the control particle, |$) c for it is easy to see that the 
change in the state of the target particle is given by 



Pout 

where 



dX P(X)U(X) Pm uHx) (3) 



4\ 



x- 2 



(4) 



with $ C (Q) = (Q|$) c is the position probability amplitude of the control particle immediately before the collision. 
In this form we see that we may regard the transformation as a mixture of classical controls with a random control 
parameter distributed according to P(X). 

As an example, consider the case for which |$ C (Q)| 2 is a Gaussian sharply peaked at a particular value Q with 
variance a. We first change the variable of integration to y = X — X, with where X = kQ. We can then factor out 
an overall unitary transformation, 



Pout = U(X) 



uHx) 



If we expand the unitary operators e ±ty ^^ h we can integrate term by term. If t&- << 1, we can truncate this expansion 
to second order. The output state in Eq.Q may be written as 

p out = U{X)p ln U\X) - I [q, [q, U(X)p m U\X)]] (5) 

o 

and r = ^fy-- We see that provided T << 1, we can treat the control as entirely classical. The correction term adds 
noise to the momentum of the target particle, as one easily sees that 

(Af) out = (Af) m + — (6) 

The mixed state in Eq.Q begs an interpretation in terms of a measurement performed on the control. Suppose 
that, immediately after the impulsive interaction between the control and the target, the position of the control is 
measured with extreme accuracy. The resulting conditional state of the target system, given a result Q from the 
measurement, is given by projecting the output state onto \Q) C and tracing out the control. Thus 

\^ (Q) )out = e- iKQ ^) in = e- iX ^) in (7) 

which indicates a classical control of the target subject to intrinsic quantum fluctuations of the control system position 
coordinate. In order to make this small we need to ensure that a is small, i.e the position coordinate of the controller 
is well defined. 

It is worth contrasting this model with how one would use this kind of impulsive interaction to make a measurement 
of the system position variable, q. We now regard the control particle not as a controller but as a measurement 
apparatus, a meter. Clearly there is no point in measuring the position of the controller immediately after the 
impulse as this does not change. Instead one must measure the canonically conjugate momentum coordinate P. For 
this to be a good measurement we need to prepare the state of the controller not with a well defined position but 
rather with a well defined momentum; exactly the opposite to what is required to use the system as a classical control. 
A quantum controller and a quantum meter are in this sense complementary devices. In the models we discuss below 
this kind of complementarity is determined by what controller variable is coupled to its environment. 
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II. QUANTUM STOCHASTIC MODEL. 

The impulsive model of the previous section is a pump to our intuition however it is rather limited. In practice 
we are more likely to be interested in a system evolving continuously in time subject to a classical controller. In 
this section we show that we can still regard the classical controller as a quantum system subject to continuous 
measurement of the controller variable coupled to the system. This of course means we need to treat the quantum 
controller as an open quantum system. It is the nature of this interaction that enables us to describe the target system 
dynamics in terms of a purely classical stochastic control. 

Consider two interacting quantum systems, a controller and a target (see figure]!]). The control system is also coupled 
to an environment in such a way that one particular degree of freedom in the environment may be monitored provding 
a continuous weak measurement of the controller variable that is coupled to the target system. The Hamiltonian 
describing the dynamics is 

H = H c + H s + HkXx + H e + H ce (8) 

where X is a (dimensionless) hermitian operator for the control and x is a (dimensionless) Hcrmitian operator on the 
system. We will assume that X has a continuous spectrum. The Hamiltonians H e and H ce represent, respectively, 
the Hamiltonian for the environment coupled to the control and the interaction Hamiltonian for the control and its 
environment. We will assume that H ce is simply linear in X so that in principle the environment can be regarded as 
a measurement apparatus for the control variable X. What we are interested in is the dynamics of the target system 
alone, given the interaction with the control system and given the entire measurement record of the control variable 
X. 

The theory of weak continuous measurement is now well establishcd|T2]- In the case of measurement of the control 
system observable X we use the theory described in |13] . There are two ways to view the measurement which we might 
loosely refer to as the inside view and the outside view. The first of these simply gives the unconditional dynamics of 
the measured system when no account is taken of the measured results. The dynamics is given in terms of a master 
equation for the measured system density operator by averaging over the state of the environment. The inside view 
gives the stochastic conditional dynamics of the measured system conditioned on a particular measurement record. The 
conditional dynamics is given by a nonlinear stochastic master equation and is sometimes called a quantum trajectory. 
The reason the dynamics is non linear is to ensure that the future measurement records are conditioned on what was 
seen in the past. 

We will take the unconditional dynamics to be given by 

^ = ~[H C + H s ,p]- in[Xx, P ] - j[X, [X, p]] (9) 

where T is the decoherence rate and determines the rate of decay of the off diagonal matrix elements of p in the basis of 
the measured quantity. It simultaneously adds noise to the observable canonically conjugate to X. The corresponding 
conditional dynamics is given by the Ito stochastic master equation, 

dp c = - l -[H c + H s ,p c ]dt - in[Xx,p c ] - [X,p c ]]dt + \J^U[X]p c dW{t) (10) 

where 

H[A]p = Ap + pA^ -tr(Ap + pA^)p (11) 

and dW(t) is the Wiener stochastic processfH] and the superscript on the state serves to remind us that this state 
is conditional on the entire measuremtent record up to time t. The classical stochastic process corresponding to the 
observed measurement record is, 

dy(t)=T{Xy + ^dW(t) (12) 

where (^) c = tr(Xp c ) is the conditional mean. The decoherence rate, T has a dual role. It determines both the 
decoherence rate and the signal to noise ratio. Indeed T >> 1 implies rapid decoherence and large signal to noise 
ratio, which we simply refer to as the good measurement limit. In what follows we set H c = so that the control 
variable has no dynamics other than what is implied by the continuous measurement made upon it. We will remove 
this assumption in the next section when we consider simultaneous measurement of canonically conjugate observables 
on the control system. 
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The conditional stochastic master equation tends to drive the system to eigenstates of the measured quantity. To 
see this let us suppose for simplicity that the coupling to the target is turned off k — > and that we ignore the free 
dynamics of the control (that is to say, we assume it is very slow compared to the measurement dynamics). The 
uncondtiional dynamics of the controller is 



and the conditional dynamics is 



d ^ = -\[X,[X, Pc ]] (13) 



dpi = ~[X, [X,p*]}dt + ^H[X]p c c dW(t) (14) 



In the good measurement limit the unconditional state will be rapidly diagonalised in the eigenstates of X. We can 
thus assume that the unconditional state for the controller, on time scales of interest, can be written as 

poo 

Pc (t) = / dXP(X,t)\X) c (X\ (15) 



Substituting this into the unconditional master equation, Eq.(|13|) we see that it implies dp c (t) — and thus P(X, t) = 
P(X, 0). On the other hand substituting this into the conditional master equation, Eq.(Jl4|, we find that the solution 
has the same form but that the conditional distribution P c (X 1 t) obeys the classical conditional dynamics 



dP c (X, t) = V2T{X - {X) C )P C {X, t)dW (16) 

In this sense the control system, in the good measurement limit, may be described as a classical stochastic process 
subject to continuous observation. Note that averaging over the noise, we have that dP(X, t) = 0, which is consistent 
with the unconditional dynamics dp c (t) — and P(x, t) — P(x, 0). This simply says that the diagonal matrix elements 
of the controller, in the basis that diagonalises the measured quantity, do not change. If we further assume that the 
initial state has a Gaussian Wigner function, the conditional dynamics will remain Gaussian and we can derive the 
following equations for the conditional mean (X) c = tr(Xpf) and conditional variance of a = tr(AX 2 pf). Taking care 
with the Ito calculus, we find that the conditional mean and conditional variance satisfy, 



d(X) c = V2TadW{t) (17) 
da = -2Tadt (18) 

The conditional mean is stochastic but the conditional variance is deterministic and decays exponentially to zero 
at the decoherence rate. In other words, for a good measurement, the conditional state will localise on a random 
eigenstate of the measured quantity. 

Turning now to the unconditional dynamics including the target system, we make the ansatz that the solution is 
well approximated, on time scales of interest to the target dynamics, by 

/oo 
dX p s (X,t)<E>\X) c (X\ (19) 
-oo 

We further assume that the controller dynamics is little changed from what it would be in the absence of the coupling 
to the target. This assumes that T is much larger than all the system frequencies. It also accords with our intuition 
for a classical controller, that is to say, a classical controller should not be affected by the system it is used to control. 
Taking the partial trace of this state over the target system gives the state of the control. Consistency then requires 
that we set 

tv sPs (X,t) = P{X,t) (20) 

We can then write the unconditonal state as 

/oo 
dX P(X t t)p.{t\X)®\X) c {X\ (21) 
-oo 

where 



(22) 
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The unconditional state of the target system alone is obtained by taking the partial trace of p(t) in Eq.(21 ) 

/oo 
dX P(X,t)p s (t\X) 
-oo 

Substituting this into the unconditional master equation we find that 



d Ps {t) = --[H s ,p s (t)]dt-iK[x,M s (t)]dt 



where 



M s (t) = 



dxXP(X,t)p s (t\X) 



To simplify this equation we use P(X, t) = P(X, 0) and approximate 

10 d 2 p(t\X) 



p s (t) « p 8 (t\X ) + 



M s (t) « X p(t\X ) + <t q 



dX 2 
3p(t\X) 



dX 



(23) 

(24) 
(25) 

(26) 
(27) 



A'o 



where Xo, 0o are the mean and variance of P(X, 0) respectively. We now proceed iteratively in powers of a. We first 
use the zeroth order solution p s (t) = p s (t\Xo) to get 



d Ps (t\X ) = --[H s ,p s (t\X )}dt-iKX [x,p s {t\X )}dt 
n 



which may be solved as 

p(t\X ) = exp 
which in turn gives the next order approximation, 



— H s t — inX a tx 
n 



Ps(0) 



-H s t + inX tx 
n 



dp s {t\X ) = --[H s ,p s (t\Xo)]dt-iKX [x,p s (t\Xo)]dt-K 2 a [x,[x,p(t\Xo) 



(28) 



(29) 



(30) 



We can pause here to note an important feature already present in the simple model of the introduction. The 
double commutator in x will cause decoherence in the position basis and drive fluctuations into the momentum of 
the controlled system. We get classical control provided we can make this term small enough over the time scales 
of interest. Provided that ctk 2 , which has units of frequency, is much smaller than the inverse time scale of the 
system, we can disregard the decoherence and take the dynamics to be very closely approximated by the Hamiltonian 
H = H s + HkXqX with Xo a classical control parameter. 

On the other hand, given that we have assumed a measurement of the controller, we do have access to the actual 
measurement record and can use this to describe the conditional dynamics of the target system. We now postulate 
that the conditional state of the total system is given by 



p c (t) = / dX P C (X, t)p c s (t\X) ® \X) C (X\ 



(31) 



where P c (X,t) satisfies the classical conditional equation, Eq.(16). It is easier to interpret this equation if we write 
P C (X, t) in functional form 



P c {X,t)= J dX t S(X - X t )T c {X t ) 



(32) 



where V c (X t ) is a distribution over stochastic paths. We can write the conditional state of the combined control and 
target system as 



Ps (t) = / dX t r c (X t ) Ps (X t ) ® \X t )c{X t \ 



(33) 
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where to ease notation we have set p s (t\X t ) = p(X t ). It is then a simple matter to see that 

dp c s (X t ) = - l -[H s ,p c s {X t )] -i K X t [x,p c s (X t )] (34) 



where X t is the observed classical control process in Eq.(16|. In the example of a Gaussian process this is given by 



the equations Eq.(18). In the case of a very good measurement X t will settle down to a constant with vanishing 
fluctuations. 

The key feature of this model is the rapid decoherence of the controller in the basis of the control variable acting on 
the system (X). This occurs because the classical controller was assumed to be continuously monitored in this basis. 
In the language of decoherence, the eigenstates of X are the pointer basis for the measurement. In general the pointer 
basis is determined by how the controller is coupled to its environment. In order to make a quantum system behave 
as a classical controller we must engineer it so that the required control variable that is to act on the system is the 
pointer variable for controller. In the simple example of this section the system-controller interaction Hamiltonian 
commutes with the pointer basis variable of the controller. 

It is interesting to contrast this result with the role of the pointer basis in a continuous measurement model[T2]. In 
the quantum theory of continuous measurement, a quantum system can play the role of the first stage of a measuring 
instrument, let us call it the meter, provided that the meter variable that is coupled to the system does not commute 
with the pointer basis variable. In the example of this section, we could use the controller as a meter by choosing 
the interaction Hamiltonian to be of the form HkPx where P is canonically conjugate to the pointer basis variable 
X. In this case a continuous measurement of X will give information on the system variable x. A quantum system 
engineered to function as a classical controller can be regarded as complementary to a quantum system engineered to 
be a measurement apparatus. In both cases a key role is played by the pointer basis of the controller/meter. In other 
words it is how a quantum system is coupled to its environment that fixes its use as either a classical controller or as 
a measurement apparatus. 

III. AN OPTOMECHANICAL EXAMPLE. 

In this section we give a less idealised example for how to engineer a classical controller as an open quantum system. 
The example is taken from the field of quantum optomechanics in which the mechanical action of the electromagnetic 
field can be used to control or measure the motion of a mechanical resonator. In particular we will consider the case 
of a single mode optical cavity field interacting with a mechanical resonator comprising one mirror of the cavity [16 . 

The interaction between the field and the mechanics is known as the radiation pressure interaction as it corresponds 
to a force applied to the mechanical element that is simply proportional to the photon number inside the cavity. The 
Hamiltonian is given by[TS] 

~2 2 

H om = hoj c a' a + 1 tt^x 2 + h(u) c l L)a) ax (35) 

2m 2 

where a, are the annihilation and creation operators for the single mode cavity field, x,p are the canonical displace- 
ment and momentum operators for the mechanical resonator and uj c is the frequency of the optical cavity mode, uj m 
is the mechanical frequency of the moving mirror with effective mass m and L is the length of the optical cavity. It 
will be more convenient to write the displacement and momentum of the mechanical element in terms of the lowering 
(b) and raising operators (tf) 



^ (i-6t) (37) 
Then, 

+ hu> m tfb + hG Q ola(b + tf) (38) 

where 



1/2 



G ° = f 9 < 39 ) 
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is known as the optomechanical vacuum coupling rate. For most systems currently in operation Go varies from a few 
Hz [17] to a fractions of a MHz [18] . We will regard the cavity field as the controller and mechanical resonator as the 
system to be controlled. 



At first sight it would appear from the Hamiltonian in Eq. ( 38 ) that the classical control variable of interest is the 
photon number inside the cavity and the required pointer variable is the photon number. However this is misleading 
as we need first to consider how the cavity is coupled to its environment. The (approximate) pointer basis variable for 
a damped cavity mode is not the intracavity photon number but rather the amplitude of the field in the cavity JL9]. 
We thus need to drive the cavity in such a way that a steady state field can build up. The small position dependent 
detuning that results from the radiation pressure interaction is then turned into a small displacement of the cavity 
field amplitude which does define the pointer basis for the damped cavity. 

The quantum theory of a single damped cavity mode is well known [T2HT3] and is described by the master equation 

[^j = -iu> c [<Ja, p] + nV[a]p (40) 

where 

V\A\p = ApA^ - - (A^Ap + pA^A) (41) 

and n is the decay rate of the photon number inside the cavity. The pointer basis states for this master equation are 
the coherent states \a) defined by a\a) = a|a)[12j[19]. In the absence of any coherent driving of the cavity, the steady 
state is simply the vacuum state |0) which is a coherent state. More interesting is the case of a cavity driven by a 
coherent laser source with carrier frequency lul- The steady state of the cavity in this case is a coherent state with 
coherent amplitude «o given by 

k/2 + iA 

where A = cu c — lo^ is the detuning of the driving laser from the cavity and units have been chosen so that \£l\ 2 are 
in units of photon flux (i.e. number per second). This result can be obtained by including in the master equation a 
hamiltonian term of the form 

H d = h*JH {a) Elc-^^ + aE* L e lulLt ) (43) 

At this point we note that we are describing the laser by a classical control parameter and this might be regarded as 
question begging. It looks like we need to go back one stage to describe the laser itself as a classical control field. This 
reasoning will lead to an infinite regress very much like the infinite regress that results if we insist on regarding every 
stage of a measurement as a quantum system. We need to put a quantum-classical cut somewhere and in the case of 
measurement this is determined by how certain key degrees of freedom of the measurement apparatus are coupled to 
the environment. However we do not face this problem here. As discussed in some detail by Noh and CarmichaelJS], 
the description of laser driving used here is good approximation given the particular kind of open irreversible system 
that is a laser. Their explanation parallels the view taken in this paper. We will not consider this further here as our 
objective is to determine under what conditions we can regard an optomechanical interaction as a classical controller 
for the mechanical resonator. 

If we now include the driving Hamiltonian and the damping of the cavity mode, the dynamics of the optomechanical 
system, in an interaction picture with respect to the cavity field at frequency is given by 

( ± = -L[H Ii p] + *D[a]p (44) 

where 

Hj = HAaU + hw m tfb + hG a^a(b + b ] ) + h{E*a + Ea r ) (45) 

and E = ^/k£l which now has units of frequency. 

We now follow the approach of the previous section; we first consider the steady state of the controller itself 
and then use this to get approximate equations of motion, including fluctuations, for the target system under the 
assumption that the irreversible dynamics of the controller is much faster than the target system dynamics. We then 
consider the conditional dynamics given that the state of the controller is continuously monitored. In the case of 
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the optomechanical model this will be done via heterodyne detection which is a means of monitoring the cavity field 
amplitude. 

The steady state of the master equation for the controller, in the absence of the coupling to the target system, is a 
coherent state given by 



Pc = |ao)(aol 



(46) 



where ao is given in Eq. ( 42 ) . In the presence of the coupling to the mechanical system, we assume that the total 
unconditional state may be written as 



P(t) 



d 2 ap b (t, a)\a) a (a\ 



(47) 



Taking the partial trace over the controller (i.e. the cavity) gives the unconditional state of the mechanical resonator 
as 



Pb(t) 



d 2 apb(t, a) 



(48) 



Taking the partial trace over the mechanical resonator gives the unconditional state of the controller as 



Pa{t) 



d 2 aP(t,a)\a) a (a\ 



(49) 



where 



P(t, a) = tr/9;,(i, a) 



(50) 



The form of Eq. ( 49 ) implies that we are defining a Glauber-Sudarshan representation for the cavity state. We may 
need to admit some unusual distributions for this function (for example, derivatives of delta-functions) in order to 
account for non-classical effects. Note that Pb(t, a) is not a normalised state of the mechanical resonator. We may 
define a normalised unconditional state for the mechanical resonator as 



Pb(t\a) = 



Pb(t,a) 
P(t,a) 



(51) 



and we have written the normalised state as if it were a conditional state conditioned on a random variable a, with 
distribution P(t, a). However this assumes that P(t, a) is a regular probability density which may not always be the 
case. 

Substituting the assumption, Eq. ( 47 1 , into the master equation, and using the correspondence rules for the P- 
representation[2Q], we find that 



dp b (t,a) 
dt 



= -iia m [tfb,pb{t,a)} - iG a \a\ 2 [b + b\ p b (t, a)} 



(52) 



+ — {(iA + k/2) + iE) p b (t,a) 



da' 



((-iA + k/2) - iE) p b {t,a) 



d d 
+iGo—a{b + tf)p b (t,a) - iG a -^a* p b (t,a)(b + b r ) 

If we take the trace of this equation over the mechanical resonator states, we find that 



8P(t, a) 
dt 



— ((iA + k/2) + iE) P(t, a) + — ((-iA + k/2) - iE) P(t, a) 
oa oa 



(53) 



+ 



d_ 
da 



(iG a(b + &t)C«)) + A (-iG a*(b + &t)(«)) 



P(t,a) 



where we have defined the conditional mean value 

(6 + 6t)(«)=tr((6 + 6t)p 6 ( t | a )) 



(54) 



The form of Eq.(53l is almost the classical Liouville for the controller density P(t,a), but the conditional mean in 
the last term gives this a more complicated dynamics, as this term needs to be found by also solving Eq.(52). The 
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classical and quantum dynamical system are tied together in a particular way. We shall return to this interpretation 
in section 



IV when we discuss the connection to work on hybrid quantum-classical dynamics. 



We can now find the master equation for the unconditional state of the mechanical resonator, the controlled system, 
by integrating both sides of Eq.(52) over the complex variable a. This takes a rather simple form after integrating by 
parts, 



= -iUmftb, p b {t)} - iG 



b + b\ I d 2 a\a\ 2 P{t,a)p b (t\a) 



(55) 



This again begs an obvious interpretation: it looks like the system is evolving with a random Hamiltonian. However 
the random variable itself is not independent as it is tied to the system through the last term in Eq. ( 53 ) . 

In order to proceed we need to make some assumptions about the rate of dissipation in the controller, i.e. the 
cavity field damping rate, n. We now assume that the cavity is strongly damped so that k >> A, Go- Under this 
assumption we can expand the total system state near the cavity vacuum state as 

p(t) = J d2a [poo(t,a)\a) a (a\ + p 01 (t,a)\a) a {a\(a - a) + p w {t, a)(a f - a*)|a) (a| (56) 

+pn(t, a) (a^ — a*)\a) a (a\(a — a)] 8^ 2 '(a — ao) 

= Poo(t, ao)|ao)a(ao| + Poi(t, a )|ao)a(ao| (a - ao) + Pw{t, a )(a ! - «o)|ao)a(ao| 
+pn(t, ao)(a f - ao)\a ) a (a \(a - a ) 

where the operators Pjk(t, a) act only on the Hilbert space of the mechanical system. This can be written in the form 
given in Eq.(47l using known correspondences [20] for the P-representation of the raising and lowering operators on 
the cavity field. We find that 



Pit) = 



d a\a) a (a\ 



(poo(t,a) + pn(t,a))(5 (2) (a - a ) 

^—p m {t,a) + ^p w (t,a) ) (5 (2) (a- a ) 
oa* oa I 



(57) 



d 2 
da* da 



Pi 



lit, a)8^ 2 \a — ao) 



The state of the mechanical resonator is then obtained by tracing out the control in Eq. ( 56 ) to get 

Pb = Poa(t,a ) + pu(t,a ) 



(58) 



If we substitute the expression in the integrand on the right hand side of Eq.(57) into Eq.(55) and integrate by parts 
we find that 



= -iu m pb, p b {t)] -iG \a \ 2 [b + tf, p oa (t, a a ) +pn{t 1 a )] 
-iG [b + tf,a poi(t,ao) + aoPio(*,a) + pn(t,a)] 



(59) 



On the other hand, substituting the approximation in Eq.(56l into the master equation, Eq.(44), and equating 
coefficients we obtain equations for These show that the off-diagonal terms rapidly approach a steady state (i.e. 

pw — Poi — 0), and we can solve for the off-diagonal terms in terms of the diagonal terms, 



Pw = 
Poi = 



iG a 
k/2 + iA 



((b + tf)p 00 -pu(b + tf)) 



iG a 



k/2 



hr(pao(b + tf)-(b + tf)p n ) 



(60) 
(61) 



This procedure is equivalent to the adiabatic elimination of the cavity field as described in |21j . Substituting these 
into Eq.(59| we find that to lowest order in 1/k, 



dpm 

dt 



[H m ,p m ]+TV[b + tf} f 



(62) 
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where 



(63) 



and H m = hu] m wb + hGo\ao\ 2 (b + b*). The target system is then seen to behave like a classically driven harmonic 
oscillator with driving determined by replacing a^a — > \cto\ 2 and a decoherence term in x. This decoherence term tends 
to diagonalise the target system in the eigenbasis of x and simultaneously adds noise to the conjugate variable p. This 
is the same result that we obtained for the simple model of the previous section. Provided that k » u m , (Goao) 2 we 
can neglect this decoherence and treat the control as entirely classical. We again see that because the control system, 
the cavity field, is an open system with a particular pointer basis (the coherent states) we can find a limit in which 
the control may be treated classically. 

To complete the discussion we turn now to the conditional dynamics that results when we continuously monitor the 
cavity field. A continuous measurement of the coherent field in a damped driven cavity requires a knowledge of two 
real numbers and can be done by heterodyne detection[12 . In essence this is a continuous simultaneous measurement 
of the two canonically conjugate quadrature phase operators of the cavity field. The resulting measurement record is 
described by two real stochastic processes, which may be combined to form a single complex stochastic process. The 
conditional master equation for a single damped cavity mode subject to heterodyne detection is 

= — iujcla^a, p]dt — i[E* a + EaJ , p c 



(dp c ) c 



]dt + nV[a]p c dt + yfim.[a]p c dZ* 



where the conditioning superoperator, H, is given in Eq.(ll), and the complex valued Wiener process dZ 
idW 2 )/V2 and 



dWi, 



dW\ 





dt 



The classical measurement record is given by the in-phase and quadrature-phase heterodyne currents, J x , J, 
obey the classical stochastic (Ito) differential equations 



(64) 
{dWi + 

(65) 
(66) 
which 

(67) 
(68) 

(69) 

Note that the superoperator, % gives zero when acting on coherent states; another indication that the coherent states 
are the pointer states for the damped cavity field. 

Given that we do have access to the measurement results we can use them to classically control the mechanical 
resonator. We now show that this is consistent with the irreversible, unconditional dynamics of the controlled system 
we have already described. With this in mind we consider the stochastic Hamiltonian 

H(t) = huj^b + hGo^^ib + tf) (70) 



dJ x (t) — n(a + a^)dt + 
dJ y {t)dt = - 



%K(a — a' 



)dt 



/2ndWi 
^V2ndW 2 



which may be more conveniently written in terms of the complex stochastic current, J — ( J x + iJ y )/2, 

dJ(t) = n(a) c (t) + ^ndZ 



a.s 



over a small time step dt the system then evolves as 

Pb (t + dt) = dUp b (t)dU' f 

where 

dU = cxp -~dHi(t) 



(71) 



(72) 



If we expand dHj(t) out to second order in the stochastic increment, taking care to keep the Ito increment dZdZ* , 
and linearise around the steady state value (a) c = cxq we find the effective stochastic increment of the unitary evolution 
operator is 



dU (t) = exp 



tfbdt 



zG |a | 2 (6 + tf)dt - i^^dWi 



(73) 



Expanding this to second order in the stochastic increment and averaging over the noise we find that the unconditional 
dynamics is given by the master equation in Eq. ( 62 ) . This indicates that it is entirely consistent to describe 



the dynamics conditionally using a classical stochastic control Hamiltonian based on the appropriate measurement 
record and the unconditional dynamics that results when we treat the open systems dynamics of the total system by 
adiabatically eliminating the controller. This is illustrated schematically in figure [l] 
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FIG. 1: Equivalent ways of describing the interaction between two quantum systems as classical control. In (A), the controller 
and the system interact unitarily but the controller is an open system subject to continuous measurement, with measurement 
record X(t). In (B) the measurement record is used as a classical control signal to directly drive the target system. In both cases 
the unconditional dynamics (average over the measurement record) and conditional dynamics (conditional on the measurement 
record) of the target system is the same. 



IV. DISCUSSION AND CONCLUSION 



We have given two examples of interacting quantum systems in which one of the systems can be regarded as a 
classical controller for the other. The conditions necessary for this are, firstly, that the controller is open in such a 
way as to define a pointer basis for the controller operator that is coupled to the controlled system and secondly that 
the open system dynamics of the controller is such that the approach to the pointer basis is very rapid compared 
to the timescales of interest for the controlled system. In the optomechanical example we saw that the key features 



are captured by two linked equations, Eq.(52| and Eq.(53|. One equation gives the unconditional dynamics of the 



controlled quantum system alone in terms of an equation for a state operator p(t, a) that depends on a classical phase 
space variable, a. The other equation gives the evolution for the classical distribution function P(a,t) that depends 
on a conditional mean of the controlled system variable occurring in the interaction (b+b^)^ a \ This structure parallels 
closely an approach to consistently combining quantum and classical dynamics developed by Diosi and others (see [7] 
and references therein). 

In [7] a classical system with canonical variables (q,p) is interacting with a quantum system. The description 
begins by defining an operator p(t, q,p) acting on the quantum system Hilbert space such that the actual state of the 
quantum system is given by 

P«C0 = / dqdpp(t,q,p) (74) 

The operator p(q,p) is called the hybrid state. The corresponding classical Liouville distribution for the classical 
system is defined by 

P c (t,q,p) = trp(t,q,p) (75) 
A hybrid observable is then defined as an operator on the quantum state space as A(q,p) with average 

(A{t)) = J dqdp tr [A(q, p)p(t, q, p)] (76) 

This is all well and good but do the dynamics correspond to anything that would make physical sense? Diosi 
shows[7] that it does provided the quantum system undergoes decoherence in just the right way and the classical 
system suffers a corresponding noise contribution. To see how this works we define a hybrid Hamiltonian H, 

H = H Q +H c (q, P ) + Hqc («, P) (77) 

where the last term is a hybrid observable intended to describe an interaction between the quantum and classical 
subsystems. The dynamics is then give by the quantum Liouville equation for the quantum system with an additional 
decoherence term. The classical evolution is given by a classical Fokker-Planck equation that includes a classical drift 
term via a noise averaged Aleksandrov bracket and an additional classical diffusion term. The quantum decoherence 
rate and the classical diffusion rate are related so as to ensure the overall dynamics is a positive stochastic process. In 
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the optomechanical example the positivity is ensured by the construction as both the classical noise and the quantum 
decoherence ultimately originate from a positive quantum evolution equation for the controller and the target system. 

The hybrid dynamics defined by Diosi and others was motivated by the problem of reconciling classical gravitation 
with quantum dynamics. We have seen in this paper that similar models arise when a quantum system is used as 
a classical controller for another target quantum system. The possibility of classical control emerges as a result of 
the interaction of the controller with its environment. This is very similar to the way in which a measurement device 
works: a classical stochastic output is possible only if the measurement apparatus is coupled to the environment in 
this right way (i.e. define a pointer basis). It would be interesting to see if a similar possibility were open for gravity. 
Could the apparent classical nature of gravity be telling us something important about the nature of the underlying 
quantum theory? Perhaps the classical gravitational degrees of freedom are coupled to hidden environmental degrees 
of freedom in just the right way to make general relativity possible. If this were true there would be an unavoidable 
decoherence in the dynamics of every quantum system moving in a gravitational field. The real question however is, 
how big is the effect? 
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